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We present a new method to derive kinetic equations for systems undergoing nonlinear transport in 
the presence of memory effects. In the framework of mesoscopic nonequihbrium thermodynamics, we 
obtain a generahzed Fokker-Planck equation incorporating memory effects through time-dependent 
coefficients. As applications, we first discuss the non-Markovian dynamics of anomalous diffusion in 
a potential, analyzing the validity of the fluctuation-dissipation theorem. In a second application, 
we propose a new ratchet mechanism in which the periodic driving acting on the particle is induced 
by the Onsager coupling of the diffusion current with an oscillating thermodynamic force. 



I. INTRODUCTION 



The characterization of systems outside equilibrium requires an a priori identification of the relevant time and length 
scale playing a major role in the dynamics. Once this task has been accomplished one then proceeds to establish the 
set of kinetic equations governing the behaviour of the probability densities and subsequently to obtain the correlation 
functions which can be contrasted with experiments. These are the main steps nonequilibrium statistical mechanical 
theories set up to analyze the dynamics of systems of very different nature evolving under out-of-equilibrium conditions 

The present-day-tendency in the study of systems working at micro- and nano-scales, based on the increasing exper- 
imental possibilities, demands the implementation of simple theoretical frameworks able to cope with the complexity 
inherent to the analysis of systems with many competing time and length scales (usually referred to as complex 
systems). Two are the main characteristics exhibited by those systems. The unavoidable presence of energetic and 
entropic barriers makes the dynamics intrinsically nonlinear _,2J- 5j. The drastic elimination of degrees of freedom 
usually performed for the sake of simplification is accompanied by the appearance of memory effects These 
two characteristics, together with the fact that frequently the system evolves under the influence of external gradients 
or fields creating an inhomogeneous environment, must be reflected in the formulation of the kinetic equations. To 
implement that task, many efforts have been devoted in the framework of kinetic theory [10|, projector operator 
techniques ^2-0 and in the theory of stochastic processes [T^-[T^. 

Two types of Fokker-Planck equations incorporating memory effects are frequently presented. The first one intro- 
duces the non-Markovianity in terms of retarded kernels, i.e., by incorporating time integrals in the equation h3|. 
In the second type, non-Markovianity is taken into account through a time dependence of the coefficients |lClj |. 

A detailed discussion explaining their differences together with experimental implications can be found 




Our purpose in this paper is to propose a new method to derive kinetic equations of the Fokker-Planck type in which 
memory effects are incorporated through time-dependent coefficients. In recent works [IJ-H^, we have shown that 
the scheme of nonequilibrium thermodynamics also applies to the mesoscopic domain in which fluctuations strongly 
influence the dynamics. The kinetic equations can be obtained from the probability density conservation law by 
specifying the phase space diffusion current which follows from the corresponding entropy production. The theory, 
called mesoscopic nonequilibrium thermodynamics, has been successfully applied to different situations involving 
transport [2^, l23| and relaxation phenomena Ell- is our objective in this paper to show its validity in systems 
exhibiting memory effects by discussing examples pertaining to the domain of nonlinear transport. 

The paper is distributed as follows. In Sec. 2, we present the derivation of the Fokker-Planck equation for 
non-Markovian systems or generalized Fokker-Planck equation. Section 3 is devoted to obtain the corresponding 
generalized Langevin equation. In Sec. 4, we apply previous results to the case of anomalous diffusion in an external 
potential. A second application to the domain of nonlinear transport in ratchet devices is carried out in Sec. 5. 
Finally, in the discussion Section we summarize our main results and comment on new perspectives. 



II. THE GENERALIZED FOKKER-PLANCK EQUATION 

Consider a system whose dynamics is determined by the set of variables a = (ai, a„); their evolution in time a{t) 
defines a stationary non-Markovian process. Let Oq denote the initial values of the variables, corresponding to t — Q. 
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The probability of finding the system in the state a at time t is then given by means of the conditional probability 
density defined as 

P(ao|a,t) = — — — — , (1) 

where P(ag;a, t) is the joint probability density and P{ai^) the stationary probability of the state ag, The 
probability density satisfies the normalization condition 



P{aa\a,t)da = l, (2) 



and is governed by the continuity equation 



d d 

-g^P{ao\a,t) -—■ Jj^aj^,a;t), (3) 

expressing probability conservation in a-space. Here Zafeoi^^i^) represents an unknown probability current. 

In order to derive the explicit ex pres sion of this current, we will follow the method proposed by mesoscopic nonequi- 
librium thermodynamics fMNET)]25l|.p^. By assuming local equilibrium in a-space, we establish the Gibbs equation 

SS{t) = -^ J fiiag,a;t)SP{aja\a,t)da, (4) 

where S{t) represents the entropy of the system, n{ag,a;,t) a nonequilibrium chemical potential and T the tem- 
perature, assumed constant. For the Markovian case, in which the joint probability factorizes in the form 
-P(Qloi2iiO = -P(£io)^(fij the Gibbs equation reduces to 

SS{t) = / f^{a,t)6P{a,t)da, (5) 



T „ 

where the chemical potential t) is related to the one of the non-Mar kovian case through 

fi{a,t) = j ^-{aQ,Q:,t)P{aQ)daQ, (6) 

and the probability f'(ag) at the initial state is normalized to the unity. 

To determine the value of the current i/c(ag,a;t) we will first calculate the entropy production cr(<). Taking the 
time derivative of Eq. Q), substituting Eq. ||2J) and integrating by parts with the condition that the flux vanishes at 
the boundaries of the system, we obtain 

If d 

o- = / i.a(ao>a;0 • -Q^lJ'{aQ,a:,t)da. (7) 

Now, following the rules of nonequilibrium thermodynamics we may establish a linear relationship between the flux 
and the thermodynamic force in the general form 

If d 

JgAao,a;t) = / Lia,a';t) ■ -^fi{aQ,a';t)da! , (8) 

where the Onsager coefficients Lij{a,a';t) are, in general, function of the state variables a and satisfy the Onsager 
reciprocal relations L_{a,a!;t) — i), where f means the transpose of a matrix 31]. These coefhcients are 

related to relaxation functions giving the correlation of the state variables with their initial conditions |30|. In the 
Markovian case, correlations decay exponentially and the L's become independent of time. In general, the L's may 
also depend on time through the conditional averages of the state variable defined by 

QL{t)= J aP{ajo\a,t)da. (9) 

To determine the expression of the chemical potential, we will use the requirement that the Gibbs equation Q has 
to be compatible with the Gibbs entropy postulate [s^l 

Sit) = -kB [ P{ao\a,t)ln^^^§^da + So. (10) 
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where ks is Boltzmann's constant, Po{a) is the probabiUty density and Sq the entropy of the reference state, considered 
at local equilibrium. By taking variations of this equation and integrating by parts, we can compare the resulting 
equation with Eq. obtaining the following expression for iJ,{aQ,a;t) 

P{an\a,t) 

fj,[aQ, a; t) = ksT in — — — h^o(a). (11) 

Here fJ.o{a) is the chemical potential at local equilibrium. If we now substitute Eq. ()ll|l into ©, we obtain the 
constitutive relation for the diffusion current in a-space 

f ~ d 
lai^Oy^yi) = - / 0.^,a']t) ■ [fesT— P(ao|a', t). 

+P{ao\a!,t)XiQL)]da', (12) 

where we have introduced the thermodynamic force 

Ma) = -^NTlnPo(a) - Mo(a)], (13) 

and defined the matrix 

1 L(a, a'; t) 

P{a.d;t)^- =Y'-' (14) 
= 1 P{aQ\a,t) 

In the Markovian limit, these coefficients must reduce to the friction coefficients in order to recover the correspon ding 
Fokker-Planck equation from Eqs. Q and (|12|) . As the L's, the /3's are invariant under time reversal symmetry 

and may also be expressed in terms of relaxation functions [SOj . Since the entropy production may in general contain 
contributions from different irreversible processes cross effects among currents and forces of identical tensorial nature 
might arise. Thus, for the sake of simplicity, the scalar product in Eq. H12I) will implicitly indicate a sum over all the 
corresponding contributions and, consequently, /3(a, a'; t) is a matrix of matrices containing the coefficients accounting 

for the possible cross effects. Now, by substituting Eq. (|12|l into (j2Jl we finally obtain the generalized Fokker-Planck 
equation 



P{ao\a,t) =0^ - J [kBTgia,a';t) ■ -^P{a^\d,t) 

+P{QL,QL\t)- X{eL)P{QLQ\QL,t)]dgl. (15) 



This equation accounts for the evolution in time of the conditional probability density for systems exhibiting non- 
Markovian dynamics when non-local effects are considered. Eq. (|15l) incorporates the memory effects through the 
time dependence of the coefficients, 'T3^,'2^]. Moreover, it generalizes the non-Markovian Fokker-Planck equation 
obtained in Refs. |0, |0 by including nonlinear contributions. 

A. The slow variable approximation 

When the characteristic relaxation time scale of the disturbances produced by thermodynamic forces is short 
enough when compared with the corresponding one in the evolution of a, the process becomes local, i.e. /3(a, a'; t) = 

(3{a;t)S{a' ~ a). Under that condition, Eq. ((T3|l reduces to 



d ^, . , d 



d ~ 

b{a;t)P{ajg\a,t) + kBT— ■ P{a;t)P{aQ\a,t) 



(16) 



which constitutes the generalized non-Markovian Fokker-Planck equation in the slow variable approximation 1] . Here 
we have defined the general thermodynamic force 6(a) by 

d ■- 

bia; t) = I3(a; t) ■ X{a) - fc^T— • /3(a; t). (17) 
= oa = 
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The Fokker- Planck equation was previously derived from the nonlinear Boltzmann equation ^3 ^tnd from the master 
equations in js^, respectively. Its Markovian version follows by considering /3(a; t) independent of time. The resulting 

nonlinear equation agrees with the corresponding one obtained in [3| from projector operators. In the case when the 
force is linear, = a, and the coefficients only functions of time, Eq. (|16|l reduces to 



d 

aP{a, t) + kBT—P{a, t) 
oa 



(18) 



where we have used the fact that P{t) — (t). In particular, for the case of Brownian motion a is simply the velocity 
of the particle. When (3ij{a;t) = f3{t)Sij, Eq. H18|l agrees with the generalized Fokkcr-Planck equation found in _15] . 



B. Determination of the coefficients /3ij 

The matrix elements Pij can be determined from the coupled set of constitutive equations (|12|l and the evolution 
equations for the relaxation functions Xij{'^i{t)) defined by 



Xu("»(0)=y aiaojp{aQ\ai,t)dai, (19) 

where p{aQ\ai,t) is the reduced probability 

p{aQ\ai,t) = J P{aQ\a,t)dai...dai^idai-^-i...dan. (20) 

The evolution equations for the Xij then derived by multiplying Eq. IjlGf) by ai and integrating over it. One 
obtains the regression laws 

f^X^J =-P^Jit)X^i&it))ao,, (21) 

For arbitrary forces, Eq. H21(l may be approximated by expanding 2LiQ.) in Q.- The coefficients are then given in 
terms of a cumulant expansion of P{aQ\ai, t). In the case of several variables, the remaining conditions for the (3ij (t)'s 
are obtained by comparison of the different terms entering the hierarchy of evolution equations for the moments of 
P(aQ|a, i), which have been previously approximated, for times t 3> I3~^{t), by the corresponding terms obtained 
from Eqs. ^ (see Refs. ^ and 

The simplest example concerns the case of a Brownian particle in which ai = u. When j3{t) — I3{t)l, with 1_ the 
unit matrix, Eq. \'2,1\ may be written as 

~Puu{t) = -Xm{i)j^Xm{t). (22) 

where we have defined Xuu{t) by using Eq. H19|l . This formula was used in [iS] to derive the linear non-Markovian 
Fokker-Planck equation. Notice that for Markovian processes x(a(i)) decays exponentially. Then, by using Eq. H21|l . 
one may verify that (3 reduces to the friction coefficient (3. 



III. THE GENERALIZED LANGEVIN EQUATION 



Our purpose in this section is to present a new derivation of the generalized Langevin equation based on our previous 
formulation. We first start from the interpretation of the probability density P(aQ|a, t) as an ensemble average of 
the density p(ao|a, t) in the a-space, i.e., P(ao|a, t) = (p(ao|a, t)). The actual value of that density differs then from 
P{aj^\a,t) in the presence of fluctuations 

6p{Q^\a, t) = p(ao|a, t) - (p(ao|a, t)) . (23) 

Consistently with this interpretation, used in the derivation of the Boltzmann-Langevin equation fsT], Eq. Q 
results after averaging the equation 
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-g^piao\a,t) = -—■^ia^,a;t), (24) 

with J.a{ao, a;t) — ^^^(ag,^;^:)^ The current $0(0.0, a; i) splits up into systematic and random contributions in 
the form 

*Jao.a;i) = (*a(«o>^;0) + *a(ao.a;0- (25) 
The random part has zero mean and satisfies the fluctuation-dissipation relation 

(ta(ao,a;t)$f (ao,a';t + r)) - 2kBTg{a,o!-T)P{aQ\a,t)6{T). (26) 

with the bracket denoting a stationary average and where we have taken into account Eq. H14() . By substituting H12() 
into (|25|l and the result into H24|l . we obtain the fluctuating Fokker-Planck equation 



d d f ~ d 



+|(a(i))X(a(t))p(ao|«',i)]da' - ^ ■ ^^ia^, a;t), (27) 

whose average coincides with the generalized Fokker-Planck equation p5|) . In the Markovian case, Eq. H27|l reduces 
to the Fokker-Planck-Langevin equation derived in [23| ■ 

The generalized Langevin equation can now be derived from the evolution equation of a{t) , obtained from Eq. Q , 
by using in Eq. (|27|l the slow variable approximation, and choosing the conditional probability p{aQ\a, t) = S{a{t) — a). 
After integrating by parts one obtains 

^^ait)^-B{a{t)) + F^{t), (28) 
where we have defined the nonlinear force term by 

B{a{t)) = P{a{t))X{a{t)) - keT-^ ■ P{a{t)), (29) 

and the random force by F^{t) — J (a, t)da. Eq. (|28|l generalizes the Langevin equation used in Refs. and 
since it incorporates the general force term (|29() . If the process is Markovian, it reproduces the corresponding Langevin 
equation since P{a,t) = /3(a) and, by using piaf)\a,t) = S{a{t) — a), one obtains /3(a(t)) — //3(a)^(a(i) — a)dcx, 0]. 
The fluctuation-dissipation theorem associated to Eq. H28|) has been derived in the linear non-Markovian case ^3 • 



{F^m^it + t)) = - (|:|(r)) {a{t)a{t + r)) , 



(30) 



whereas in the nonlinear Markovian case one has 35] 

{F^{t)F"it + r)) = {aB{a) + B{a) a) 5[t), (31) 

The last equation has been obtained for an arbitrary dependency of the stochastic force on a. 

The case analyzed in i.e. XJ^qit)) — a{t), follows easily from our formalism. From Eqs. H28|) and H29|l one 
obtains 

j^a{t) = -^{t)-a{t)+F^it), (32) 

or equivalently to the generalized Langevin equation [Tll | 

d 



a(t) = - I3{t~ s) ■ a(s)ds + F"'(t), (33) 
at Jo = 

where f3{t) is the memory function. To prove that equivalence, and following |0| we will first take the average of Eqs. 
(|32|l and (|33|l . By using Laplace transforms, the solution of Eq. H33|l is given by a{t) = x(^)£^0' where the relaxation 

function is x(f) — L^^ |[s + /3(t)]^^|. Taking now the time derivative of a(t), eliminating and using the definition 
, we finally arrive at Eq. H32|l . 
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IV. BROWNIAN MOTION IN AN ANHARMONIC HEAT BATH 

In this section, we will analyze the nonlinear dynamics of a Brownian degree of freedom in an anharmonic bath. 
The state variable will be in this case a = (r, u), where r is a coordinate and u its conjugate velocity. 
Let us consider a reference local equilibrium state characterized by the stationary probability density 

where /ig is a reference chemical potential per unit mass and (j){'r) the potential per unit mass imposed on the system 
by the heat bath. For this system, the corresponding Gibbs equation is 

11 If 

Ss = -5e + -pSp^^ ^ T ("^^^ 

Here s and e are the entropy and the internal energy per unit mass, respectively, p the hydrostatic pressure, cs = 
the Brownian mass fraction and p = p{uQ,ro\r,t) is defined through Eq. 1)20(1 . By taking the total time derivative, 
= ^ + ■ V, of Eq. P5|l and using Eqs. (|M)l and (jSH), we arrive at the entropy production 



du. (36) 
The average velocity v is defined through the first moment of the probability density 

pv{'r,t) = m / uP{uo,'ro\u,f,t)du. (37) 



The entropy production a adopts the usual form in terms of flux- force pairs. The two first contributions result from 
the diffusion current in -u-space and the current J = {u — v)P coming up when inertial effects are relevant. The 
average contributions arising from the pair J and V(/)(r) vanish by virtue of Eq. 1(37(1 . j3Cl| |. The third contribution 
is genuine of the effect of the potential and introduces the current = uP. In the Markovian case, when external 
forces vanish, the second term on the right-hand-side of l(3t)() contributes to the entropy flux, reproducing in this form 
the results of Ref. (2^. The constitutive laws giving the currents J^, J and in terms of the forces are, respectively 

Ju = -hu{t) \uP+^—\- Pur{t)^VP - h^{t)PV<j), (38) 



/ = -Pruit) Up+^^]- hr{t)^VP - Pr4t)PVct>, (39) 

m ou ) m 

and 

= ~ku{t) i^P +^^}~ ^^r{t)^VP - Mt)PVcb, (40) 
1^ m ou ) m 

where the Onsager coefficients have been defined consistently with Eq. ((14(1 . By substituting Eq. ((38(1 into the 
corresponding continuity equation, one obtains the generalized Fokker-Planck equation describing the non-Markovian 
dynamics of the Brownian degree of freedom in the presence of arbitrary bounding forces 



dP d 
— + V • = — • <; I3auit) 
at ou ' 



uP 



keTdP 
m du 



- Pnc^{t)X{r)P + M^/3sH0Vf1 , (41) 
m J 

where X{r) ~ — V0. In the linear case, it reduces to the corresponding one obtained through the conditional 
probability associated to the Langevin equation used in M 

As a consequence of the existence of local equilibrium 
in u-space, f3su{t) multiplies both the drift and the diffusion terms entering ((41(1 . The presence of forces adds two 
terms that make the fluctuation-dissipation theorem not longer valid. For Markovian processes in which (3su{t) = /3 
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, Pucp{t) = 1 and ^ur(t) — Oi Elq. l|lT)l reproduces the corresponding Fokker-Planck equation 0. By integrating Eqs. 
and H4Q(I over u and combining the resulting equations, we obtain the diffusion current Joif^t) = pv(f,t) 



Jd 



Vp~pX{r} 



(42) 



where the existence of local equilibrium imposes fi{t) ee /3^^[f3^^{t) - (3p^{t)(3p}{t)] = f3^^[f3^p{t) - f3pp{t)(3p}{t)]. By 
averaging Eq. H41|l over u-space and substituting Eq. (|42|l into the resulting equation, we obtain the non-Markovian 
Smoluchowski equation for the density field p 



77(t)V • [pXif) j 

where we have defined the diffusion coefficient D{t) in terms of the generalized Einstein relation 



D{t) 



(43) 



(44) 



from which it follows that the slow modes of non-Markovian systems satisfy a modified fluctuation-dissipation relation. 
For the particular case of a harmonic potential, (j){r) = ^w^r^, where w is the characteristic frequency, Eqs. (|4f |l and 
()43|l reproduce the generalized Fokker-Planck and Smoluchowski equations derived in and 0|, respectively. Now, 
multiplying Eq. I|4f(l by u and integrating over it, we obtain the evolution equation of the average velocity v 



Pj^ = -Puu[t) 



pv - f){t)pX{r) + D{t)Vp 



(45) 



where we have assumed that the pressure tensor of the ideal Brownian gas is given by P = pi. For times satisfying 
T 3> we have imposed that Eq. must coincide with Eq. obtaining the conditions I3^^{t) = I3^u{i)f}{t) 

and (3urit) = I3u<i){t) ^ Ij consistent with the corresponding ones for the Markovian case. /3sa(i) and f]{t) may be 
obtained by using the regression law lEU and Eqs. EH) and 

When local equilibrium in r-space is not fulfilled, the coefhcients appearing in the diffusion and drift terms of 
Eqs. (|39|l . H4U|) and (|42|) must be found independently. Then, even in the absence of external forces, the fluctuation- 
dissipation theorem (|44|l is no longer valid. This fact occurs in systems exhibiting anomalous diffusion, which are 
described by non-Markovian Fokker-Planck and Langevin equations, and 0|. 



V. RATCHET EFFECT INDUCED BY ONSAGER COUPLINGS 



In this section, we will apply our formalism to show a new manifestation of the ratchet effect which originates from 
the Onsager coupling of the probability current and a periodic thermodynamic force, when fluctuations are rectified 
by an external asymmetric potential. We will show that in the scale considered the non-Markovian dynamics makes 
the ratchet effect more pronounced. 

The continuity equation in the presence of external forces takes the form 

dP ^ dP d ^ 

^ + V-^P + F~ = -^. J„^, (46) 
at ou ou 

where the external force per unit mass is F{r) = — Vy(r), where V{r) represents the ratchet potential. In the present 
case, we will consider that the probability density of the reference state is given by Eq. H34(l with (/)(r) = 0. 

Inhomogeneities in the temperature of the bath induce a heat flux Jq whose appearance modifies the entropy 
production of the system. In this case, the Gibbs entropy postulate (|f 0|l and the Gibbs equation H35|l must be 
complemented with the balance equation for the internal energy e 



De 

In 



-V-J„ (47) 



where viscous effects have been neglected. By considering the Gibbs-Duhem relation |30j , the entropy production can 
be written in the form 

a=-j;.i,Vr-^/j..|dz/, (48) 
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where the modified heat flow J* is given by 



J: = J a 



1 



m I u J du. 



(49) 



The second term on the right-hand-side of Eq. H49|l represents a local flow of kinetic energy. 

The diffusion current Js is, according with Onsager theory coupled to the temperature gradient and is given by 



keTdP 
m du 



PuT(t)P 



VT 



Similarly, the heat current can be formulated as 

j; = -pTT{t)VT ~ pTu{t)pv, 



(50) 



(51) 



where f3TT{t) = '^^a^*'' , and the coupling coefficient Pxa has been defined accordingly with Eq. (|17|) . By substituting 
Eq. H50|l into the corresponding continuity equation for the probability density, we obtain the generalized Fokker- 
Planck equation 



— + V-UP + F- — 



dP 

'di 
d_ 

du 



dP 



uP 



du 
ksTdP 



du 



+ PuT{t)P 



(52) 



which describes the transport of a Brownian particle in an thermally inhomogeneous environment. 

By multiplying Eq. H52|) by u and integrating over it, we obtain the corresponding evolution equation for the 
momentum 



= -Puu{t) [pv + D{t)\Ip +DT{t)pVT] + pF{r), 



(53) 



where Puu{t) has been given in Eq. I|22|l and we have defined the diffusion coefficient D{t) — !^^/3^i(t) and the 

thermal diffusion coefficient Drit) — ^1 + jT^Paritij ■ For convenience, we will assume that D{t) and Drit) 

will be only functions of time. The diffusion current Jd is obtained from Eq. (|53|l at times t ^ Pm^t)- By following 
the steps indicated in Sec. 4, we arrive at the generalized Smoluchowski equation 



1^ = D{t)v^p + Drm ■ ip^T) - ■ [pm 



(54) 



where memory effects enter the equation through the time dependence of D(t), Drit) and Puu{t). Since in practice the 
temperature inhomogcneities originated by density gradients are negligible j33, by substituting the diffusion current, 
implicit in Eq. (|54|l . and Eqs. (|49|) and H51|) into H47|) . we arrive at the heat conduction equation 



dT 

'~dt 



(55) 



where c„ is the heat capacity. We have used the relation, e = c^T and defined the effective heat conduction coefficient 



Xeffit) = (iTT{t) l + pDT{t) 



Pj!T(t) 
l3TT{t) 



To obtain Eq. we have also neglected contributions arising from the 



kinetic energy fiow |23|, and assumed temperature gradients such that ( )" <C 1, where Iq is a characteristic 

length of the bath. For I3tt{^) large enough (~ IQ'^gr cm s^^K'^) and 13^(0), /^st(0) ^ 1, Xeff is practically 
constant in first approximation. 

To illustrate the ratchet effect, we consider the case of a Brownian particle moving in a bath, filling the space 
X > 0, in which a periodic temperature gradient is established when T(0, t) = To[l + cos{uot)]. This situation may be 
accomplished when the distance of the particle to the wall is smaller than the wave length of the thermal waves. In 
such a case the spatial dependence in the propagating wave 



T{x, t) = To [1 -I- e^*''-' cos{kx - cot)] , 



(56) 



solution of Eq. H55() . can be neglected. Here w is a frequency and k = (^x^)^ ^ wave number. The temperature 
gradient is then given by VT = kTQ[sm(ujt) — cos{Lot)]. This approximation is valid for values of the parameters 
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indicated in the figure captions, all of them accessible in experiments. The particle in turn is under the influence of 
the asymmetric potential (see Fig. 1). We will assume that it undergoes anomaluos diffusion with D{t) ^ and 
Puu{t) t~2 . This particular dependence has recently been observed for the Brownian motion of a particle inside an 
eukaryotic cell (||. The diffusion coefficient D{t) has been obtained through the time derivative of the mean square 
displacement measured in 0, and (3ss{t) by using Eq. (|^ . The thermal diffusion coefffcient Drit) is then determined 
by assuming that the Soret coefficient St = "tj^ is constant in first approximation. According with experimental 

results [s^, we will consider that -0(0) ~ lO^^cm^s^^ and that St lO^^A'"^. Then, by using p{x, t) = S{x{t) — x) 
in Eq. H53|l . after integrating over the system volume we obtain the evolution equation for the position of the particle 

+ F'^{t). (57) 

x—x(t) 

The inertia term is maintained since non-Markovian effects arise at short time scales. 

Positions and velocities of the Brownian particle have been determined from this equation by using a fourth order 
Runge-Kutta method. In Fig. 2, we show the displacement as a function of time for a given initial condition and a time 
periodic force. After a transient regime determined by (3fi^(t), in the non-Markovian case the current, proportional 
to the slope of the curve, becomes preactically constant. The average motion has the same frequency as the time 
periodic force. 

In Fig. 3, we represent average displacements for different values of the parameters. In curves a) and b) the strength 
of the potential and the frequency of temperature gradient have been used to control the magnitude of the current. 
We have verified that no net current appears when the potential is symmetric. As a consequence of the fact that the 
damping coefficient decreases with time, we may conclude that, under similar conditions, the non-Markovian system 
is more efficacious in rectifying the motion than the Markovian one. 

VI. DISCUSSION 

In this paper, we have proposed a new method to derive the generalized Fokker-Planck equation describing the 
dynamics of systems when memory effects are relevant. Our formalism is able to derive kinetic equations from the 
probability density conservation law by providing expressions for the probability current obtained from the entropy 
production in phase space. It permits to incorporate in a simple way basic ingredients proper of systems outside 
equilibrium as the presence of external forces and of inhomogeneities in the bath. 

We have presented two applications. In the first one, we have discussed the non-Markovian Brownian motion in a 
bath exerting a force on the Brownian particle 14J, 15] . The Fokker-Planck and Smoluchowski equations have been 
derived for an unspecified force by avoiding the difficulties inherent to the habitual method, obtaining them from 
the corresponding Langevin equation. Our method relates the validity of the fluctuation-dissipation theorem to the 
existence of local equilibrium in the space of the mesoscopic variables. An example in which this feature becomes 
manifest is the case of anomalous diffusion 0,0 ■ In the second example, we have analyzed the role played by the 
presence of inhomogeneities in the bath, due to an imposed temperature gradient, in nonlinear transport. We have 
found a new ratchet effect which originates from the cooperation of a thermodynamic periodic force with the nonlinear 
potential. The positiveness of the entropy production imposes restrictions on the values of the parameters entering 
mathematical models. 

These results confirm the usefulness of our method which offers a very promising framework in the study of systems 
operating at the mesoscale for which the presence of nonlinearities, the participation of many degrees of freedom and 
the inffuence of inhomogeneities in the environment become a common feature of the dynamics. 
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FIG. 2: Position of the Brownian particle versus time. Dotted lines represent different realizations of the noise, whereas the solid 
line corresponds to the average displacement of the Brownian particle for 20 realizations of the noise in the non-Markovian 

case. The dashed line represents the average position in the Markovian case. The values of the different parameters are: 
F(0) = 1 X 10"^firr cm^s"^ VT ~ 1.3Kcm'\ w = 9 x 10"'^s"\ fc ~ 2.1 x lO^^cm^^ and m ~ IQ-^gr. 
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FIG. 3; Average displacements for different initial positions, frequencies and strengths of the temperature gradient with -F(O) = 
6.8 X lQ~^gr cw?s~'^, and m ~ 10"''gr. In cases a) and c): VT ~ 1.3Kcm~^ and lo = 9 x 10~^s~^and A: ~ 2.1 x 10"^cm"^ 
In cases 6) and d): VT ~ 0.2Kcm~^, u; = 6 x W^s'^and fc ~ .1 x 10" cm . Curves a and 6 correspond to non-Markovian 
cases, whereas c and d result from Markovian dynamics. 



